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Abstract. We present an instability occurring in the in- 
ner part of disks threaded by a moderately strong vertical 
(poloidal) magnetic held. Its mechanism is such that a spi- 
ral density wave in the disk, driven by magnetic stresses 
(rather than self-gravity as in galactic spirals), becomes 
unstable by exchanging angular momentum with a Rossby 
vortex it generates at its corotation radius. This angular 
momentum can then "leak" as Alfven waves emitted to- 
ward the corona of the disk providing, as an element of the 
accretion process, an energetic source for a wind or a jet. 
As galactic spirals, this instability forms low azimuthal 
wavenumber, standing spiral patterns which might pro- 
vide an explanation for low-frequency QPOs in low-mass 
X-ray binaries. 
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1. Introduction 

The accretion disks of objects such as Active Galactic Nu- 
clei, X-ray binaries or Young Stellar Objects are very com- 
monly observed to emit jets, whose formation seems to be 
inherent to the accretion process. For instance in YSOs 
the power of the outflows has been found to be correlated 
with the accretion flux (Cabrit and Andre, 1991 ; Edwards 
et al. , 1993). In the micro-quasar GRS 1915 recent multi- 
wavelength observations (Mirabel et al. , 1998) show the 
ejection of plasmoids coinciding with the infall of the inner 
region of the disk toward the central black hole. It is thus 
considered as a strong point in favor of MHD models of 
jets that, from their very beginning, these models found 
that the jets are very efficient at carrying away the angular 
momentum extracted from the disk, and thus at permit- 
ting accretion to proceed. This has led to the success of 
a family of accretion-ejection models for these objects ( 
see e.g. Blandford and Paine, 1982; Pelletier and Pudritz, 



1992). 
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On the other hand existing models of turbulent accre- 
tion in disks always involve a transfer of angular momen- 
tum which is radial, i. e. within the disk: this is built into 
the ad- hoc model of a disks (Shakura and Sunyaev, 1973), 
where turbulence results in a viscous radial diffusion of the 
angular momentum. This is also true of specific models of 
disk instabilities (Papaloizou and Pringle, 1984; Tagger 
et al. , 1990, hereafter paper I; Balbus and Hawley, 1991) 
which all rely, in one form or another, on a radial exchange 
of angular momentum. This may appear as a contradic- 
tion between disk and jet models. We believe that it is one 
of the main reasons why it has always been very difficult, 
though not totally impossible, to develop models connect- 
ing jet solutions to disk models (Ferreira and Pelletier, 
1993; Ogilvie and Livio, 1998), in such a manner that the 
angular momentum of accretion is redirected toward the 
jet. 

We present in this paper a new instability mechanism 
which might help solve this contradiction. The instabil- 
ity appears in the innermost region of the disk (a few 
times the inner radius), where roughly half of the accre- 
tion energy and angular momentum must be disposed of. 
Its physics is such that the angular momentum extracted 
from the disk is not emitted radially but will ultimately 
end up as an Alfven wave traveling to the disk corona, 
where it can thus power the jet. 

This results from the combination of three ingredients: 

— We have shown in previous work (Paper I) that spi- 
ral waves driven by magnetic stresses (rather than 
by self-gravity as in galactic disks) propagate in disks 
threaded by a vertical (poloidal) magnetic field. One 
can thus consider that the family of spiral den- 
sity waves in disks concerns disks with an attractive 
long-range force (self-gravity) as in spiral galaxies, 
disks with only the local pressure force (giving the 
Papaloizou-Pringlc instability), and disks with a repul- 
sive long-range force, resulting from magnetic stresses. 
These waves can be amplified in the region of their 
corotation radius. The amplification is strong for self- 
gravity driven waves, and very weak in the other cases. 
However the Papaloizou-Pringle instability concerns 
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only waves with very small wavelengths (of the order 
of the disk thickness), whereas self-gravity or magnetic 
stresses move the unstable range to larger scales (small 
azimuthal wavenumber). 

— Our result, as well as most of the theoretical work on 
the propagation and amplification of spirals, was ob- 
tained in the classical shearing sheet model, which al- 
lows a thorough analysis of the physical processes in- 
volved but neglects an important effect: the gradient 
of vorticity in the equilibrium flow of the disk, and as 
a consequence the corotation resonance, whereby the 
density waves couple to a vortex at their corotation 
radius []; in a paper to be published independently we 
will present an extension of the shearing sheet model, 
keeping into account the vorticity gradient in the disk, 
and thus the physics of the corotation resonance. It 
will allow us to analyze in more details the coupling to 
this vortex, which is in fact a Rossby wave, associated 
with the vorticity gradient in the disk; the propaga- 
tion of the Rossby wave is quenched by differential ro- 
tation, but this still allows it to exchange energy and 
momentum with the density wave. In the present pa- 
per we will only use the exact cylindrical geometry, so 
that the physics of Rossby waves is fully taken into 
account. 

The effect of this resonance is only minor for self- 
gravity driven spirals, which are very unstable anyway 
(Pannatoni, 1983). It does not change much the order 
of magnitude of the growth rate for the Papaloizou- 
Pringle instability (Papaloizou and Pringle, 1985; 
Narayan et al. , 1987). It has not been considered thus 
far for the magnetically driven spirals. 

— On the other hand, we should expect a strong influ- 
ence in the magnetically driven case: in the Papaloizou- 
Pringle instability, density waves are evanescent in the 
corotation region and are thus exponentially small at 
the corotation radius, where they couple to the vortex. 
In the magnetically driven case, density waves create 
magnetic perturbations which decrease only weakly in 
the corotation region. This allows their coupling to the 
vortex to be much more efficient. More precisely we 
will find the coupling algebraically, rather than expo- 
nentially, small, the small parameter being essentially 
m (the azimuthal wavenumber) times the ratio of the 
disk thickness to radius. We will find this sufficient to 
get a sizable amplification of the density wave. 

We will show that this instability forms normal modes, i. e. 
standing, exponentially growing patterns. They are local- 
ized close to the inner radius of the disk, and thus cannot 
be considered to cause accretion in a large portion of it. 
However in this inner region, where roughly half of the 

1 the radius where the angular phase velocity of the wave 
equals the angular velocity of the gas, not to be confused with 
the radius where the gas rotates at the angular velocity of the 
central object. 



accretion energy and angular momentum are generated, 
they can cause accretion and energize a jet. 

One can have two views of this amplification process: 
we can consider it, as we did above, as the amplification 
of the spiral wave by interaction with the Rossby vortex 
it generates. But we can also concentrate on the evolution 
of the vortex: differential rotation leads it to be sheared 
away as time evolves, so that an initially leading pattern 
becomes strongly trailing; but at the same time the vortex 
generates a spiral density wave which, after reflection at 
the inner radius of the disk, returns to regenerate it as a 
leading feature, starting a new cycle. 

We present two complementary computations of the 
amplification due to this resonance. The first approach 
is analytical, in the form of a variational principle which 
allows us to discuss the physics, the sign and the order 
of magnitude of the amplification. We find in particular 
that the parameter controlling the sign of the resonant 
effect is not as usual the gradient of the specific vorticity 
in the equilibrium flow, but a different quantity involving 
the magnetic field strength. This first approach allows us 
to extract and describe in details the physical processes 
involved in the formation of normal modes and in their 
amplification. 

In a second step we present numerical solutions which 
confirm this behavior, and allow us to exhibit exponen- 
tially growing normal modes affecting the inner region of 
the disk. Their growth time is typically a fraction of the 
rotation time of the disk at the corotation radius. 

We derive both of these results in the simple model of 
an infinitely thin disk in vacuum. In this case the energy 
and angular momentum extracted from the inner region 
have to stay in the disk: they are stored in the vortex at 
corotation, rather than carried away by a density wave 
travelling outward as in the usual amplification process 
of spiral density waves. In the Papaloizou-Pringle case, 
Narayan et al. (f987) have argued that the accumulation 
of angular momentum in the vortex might eventually sat- 
urate at a finite amplitude, thus ending the amplification 
process. We will briefly discuss how here, if the disk were 
more realistically embedded in a low density corona, the 
angular momentum would be emitted as an Alfven wave 
emitted vertically along the field lines. We will reserve 
the detailed discussion and computation of this effect to a 
forthcoming publication. In our conclusion we will discuss 
possible astrophysical consequences of this new mecha- 
nism. 

Throughout this paper we will consider normal modes, 
i.e. exponentially growing perturbations characterized by 
a single frequency lo. As usual in disks with differential ro- 
tation, the amplification mechanism would also apply to 
transient shearing perturbations, which become normal 
modes only when a reflexion at the inner radius allows 
them (if they obey an integral phase condition) to start 
a cycle of wave reflection and amplification. The ampli- 
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fication process is the same. However it is weak enough 
that the transient amplification of perturbations starting 
at the noise level is probably insufficient to cause accre- 
tion; thus only normal modes, whose exponential growth 
allows them to reach large amplitudes, should really mat- 
ter in accretion disks. 

Curry and Pudritz (1996) have performed a modal anal- 
ysis of an MHD disk, in an incompressible MHD cylinder 
limit. Incompressibility did not allow them to make a con- 
nection with the magnetized spirals of paper I, and thus to 
identify the effects discussed here. They study in details 
the Alfven resonance ( u — mfl(r) = k z VA, where k z is 
a vertical wavenumbcr), which is very close to corotation 
where our main effect occurs (corotation is not a singu- 
larity in their approximation). It is very likely that the 
Alfven resonance is a form of the coupling to vertically 
propagating Alfven waves which will be briefly discussed 
below. However in the context of our modes, which are 
essentially constant over a disk scale height, k z looses its 
meaning and this would need to be studied with a more 
elaborate treatment of the vertical dimension, giving the 
vertical ejection of angular momentum associated with our 
modes. We defer this to a future publication, and will 
present in Appendix C only a brief sketch of the coupling 
with Alfven waves far from this resonance. 

2. The basic model 

2.1. Thin disk geometry and relation with other 
instabilities 

We restrict ourselves to perturbations of an infinitely thin 
disk, neglecting vertical velocities and variations of the 
perturbations across the disk. This means that we will 
not discuss the magneto-rotational instability (Velikhov, 
1959; Chandrasekhar, 1960; Balbus and Hawley, 1991), 
which relies on the existence of a wavelength in the ver- 
tical direction. In a previous work (Tagger et al. , 1992) 
we discussed the vertical structure of disk perturbations, 
and showed that there always exists a solution with weak 
variations across the disk; this solution reduces, to lowest 
order in the disk aspect ratio (thickness/radius), to the 
infinitely thin solution considered here. 

Thus we will totally discard here the vertical struc- 
ture of the wave within the disk, which is essential to the 
magneto-rotational instability. There are two reasons for 
this: one is that these instabilities are entirely distinct 
and proceed from very different mechanisms; the second 
reason is that the magneto-rotational instability occurs 
only when the magnetic pressure is weak, i.e. the ratio 
(i = 8np/B 2 of the thermal to magnetic energy density is 
larger than one whereas, as will be seen below, our insta- 
bility is significant only when f3 becomes of the order of 
one, or weaker. We will thus consider disks where the mag- 
netic pressure is strong, in the sense that (5 1, but still 
much weaker than the gravitational energy so that mag- 
netic support of the disk remains weak. When magnetic 



support becomes strong our instability becomes the inter- 
change mode studied by Spruit and coworkers (Spruit and 
Taam, 1990; Spruit et al. , 1995) in disks with and without 
differential rotation. Indeed our basic assumptions are es- 
sentially the same as theirs. In an approximate analysis of 
the effect of differential rotation Spruit et al. (1995) found 
that it makes the amplification of the interchange mode 
only transient, and negligible unless the magnetic field is 
very strong. In fact the vortex part of our instability may 
be viewed in some respects as an interchange, since the 
vortex effectively exchanges flux tubes across the coro- 
tation radius, and the amplification relies on the radial 
gradients of Bo and S (the surface density) - although in 
a different form than that found by Spruit et al. . The cou- 
pling between the vortex and density waves, which allows 
the formation of exponentially growing normal modes, al- 
lows our instability to act already at (3 <~ 1, i.e. at much 
lower magnetic pressure than in their case. A more de- 
tailed comparison is difficult, because we use a slightly 
different model: Spruit et al. consider magnetic fields sup- 
ported by equilibrium currents only in the disk, resulting 
in a jump of B r across the disk, whereas for simplicity we 
assume that the magnetic field gradient at equilibrium is 
only due to external currents (e.g. ring currents inside the 
inner disk radius), so that the equilibrium field is purely 
vertical across the disk. 

Our justification for this is that we are interested in a 
regime where the interchange instability should be negli- 
gible. On the other hand in more recent work, Stehle and 
Spruit (1998) found in non-linear simulations strong in- 
stabilities even in situations which should be stable to the 
interchange mode. It is possible that this is a manifesta- 
tion of the instability we present here (in particular the 
density and magnetic field profiles they use should give in- 
stability, according to our criterion). But the growth rates 
they find are rather large, and the modes have a more 
global structure than we would expect. Although no def- 
inite conclusion can be reached at this stage, we rather 
believe that their result is due to the reflective boundary 
condition they use at large r: in the Papaloizou-Pringlc 
case (Narayan et al. , 1987) and presumably also in the 
magnetic case, this leads to a strong form of the spiral 
instability, due to the formation of a double cavity be- 
tween the corotation radius and the inner and outer radii 
of the disk. This boundary condition, often used in disk 
theories and simulations, can thus be quite misleading, al- 
though one might look for astrophysical justifications in 
specific contexts, if some mechanism can act to truncate 
the disk at large radii. Otherwise a radiation condition at 
the outer radius must be used, as we do here. Simulations 
with different boundary conditions, or increased outer ra- 
dius, would then be necessary to fully understand their 
result. 

Our simple model of a disk in vacuum also means that 
we cannot fully justify, in the limits of the present paper, 
our claim that our instability may eventually redirect to- 
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ward the disk corona the angular momentum it extracts 
from the disk. This can however be easily understood if 
one remembers that the disk is threaded by magnetic field 
lines, along which Alfven waves can propagate as soon as 
there is any small plasma density in the corona of the disk; 
the vortex at the corotation radius moves the footpoints 
of the magnetic field lines, and this motion (of torsional 
nature, whereas the motion associated with the spiral is 
essentially compressional) will act as an initial condition 
for the emission of Alfven waves. We can thus consider 
that, although this is not contained in the present descrip- 
tion where the disk is in vacuum, the energy and angular 
momentum deposited in the vortex by the instability will 
end up in these Alfven waves propagating in the corona. In 
Appendix C we show how this occurs in the region where 
density waves propagate, in the WKB approximation. The 
full computation, valid in the corotation region, is much 
more complex and we defer it to a future publication. 

2. 2. The magnetic field 

We consider a disk threaded by a vertical magnetic field. 
For simplicity we discard radial or azimuthal (toroidal) 
components of the equilibrium magnetic field, and assign 
its radial dependence to external currents^]. Thus at equi- 
librium the disk is characterized by its surface density S(r) 
and magnetic field B = B$(r)e z . 

The disk lies in vacuum: then the absence of currents 
allows us to describe the perturbed magnetic field outside 
the disk by a magnetic potential: 



(2) 



B = -sign(z)V$ A/ 
with 

A$ M = 0. 

We work in cylindrical coordinates (r, <&, z). The sign of z 
in Eq. (Q) ensures that, if <E>m is even in z, the perturbed 
vertical field is also even, while its horizontal (r and #) 
components are odd. From their change across the disk 
we get the currents: 

jr = 2^ M S(z) 
d 

jo = -2—§ M 8{z) 

and the horizontal magnetic stresses acting on the disk: 



n 



w -2B ^ M 8(z) 
or 

d 



Fjj = -2B ^ M 6(z ) 

2 This choice is made in order to trim down our model to 
the minimal necessary physics, for the sake of clarity in this 
paper where we wish to present and discuss the basic physics 
of our instability. Taking into account equilibrium currents in 
the disk, which cause a finite inclination of the field lines at its 
surface, is straightforward, and preliminary results show that 
they increase significantly the growth rate of the instability. 



On the other hand the vertical component of the field: 
d 

B z = -sign(z) — $a/ 
oz 

gives at the disk: 
oz z 

where B® is the perturbed field in the disk. Thus through- 
out space we can write: 



A$ M = -2B°8{z) 



(3) 



This can be compared with the Poisson equation giving 
the perturbed gravitational potential in a self-gravitating 
disk: 



A$ = AnGa8(z) 



(4) 



where a is the perturbed surface density. We can thus 
compute the magnetic potential by the classical Poisson 
Kernel (Binney and Tremaine, 1987): assuming perturba- 
tions varying as exp(imd) we have: 



*M(r) = / dr'K m {^)B°{r') 



where 



K m (x) 



dd cos 



tt J Q (x 2 - 2a; cos + l) 1 / 2 



(5) 



(6) 



n\ 2.3. Perturbation equations 



We write the linearized Euler equations for perturbations 
varying as e l { m ®- ut ) : 



2flv^ 



El 
s 



, TT/ F$ 

— ILUVfl + VVV r = 

Xj 

where the forces in the right-hand sides comprise pressure 
and magnetic stresses. Here Cj = to — mfl(r), £1 is the rota- 
tion frequency of the equilibrium flow, and W is known in 
the context of galactic disks as Oort's second constanifl: 



W = 



2Q 



1 d 

r dr 



(r 2 n) 



K is the epicyclic frequency, and it is worth noting already 
at this point that W is the vertical component of the vor- 
ticity of the equilibrium flow. 

We find it more simple to change variables to: 



s = bar 

U = rv r 
V = rv§ 



3 we will not use here the more classical notation B, in order 
to avoid confusion with the magnetic field 
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We get: 



where 



iu>U - 2QV 



2 d , B d ^ 



'ds 



AujV + WU = -imc 2 s h - 2im-£-$ M 



where we have defined 



E ds 
Bo 
E 



s' 



(7) 
(8) 



T ~ 2 2 

L = uj — k 

Substituting into Eq. (|l(i|), and using: 

o 

W = 2Q + Q', —Co = -mfi' 
os 

where the prime denotes the derivative with respect to s, 
we get: 



a is the perturbed density, and cs is the sound speed. 
For simplicity we have assumed an isothermal equation of 
state, and we will consider throughout this paper that cs 
is constant. We also write the continuity equation: 



i) 



iujr 2 a = ™(EU) 
os 



imSV 



(9) 



Our system of equations is closed by the vertical compo- 
nent of the induction equation : 



iujr 2 B° = -JLfBoU) - imB V 
os 



(10) 



which can be usefully compared with Eq. (j^) : they corre- 
spond respectively to the conservation of magnetic flux 
and total density in a flux tube. Let us now consider 
Eqs. |||l0|) and (@-|). The parallel between a and B z , 
$ and &m, shows that, as explained in paper I, the mag- 
netic field appears in the mathematical form of a negative 
self-gravity, i.e. & long-range force whose action is repul- 
sive (the magnetic field stiffening the field lines and thus 
resisting, rather than favoring, the formation of overden- 
sities). We will hereafter, for simplicity, assume that the 
self-gravity of the disk is negligible but its inclusion in 
our description to discuss more massive disks would be 
straightforward. 



il>r 2 

os EL os 



A I 



d ( MmVLB 2 



ds V EL 
After some algebra we write this as: 



_ 9 B 2 uj 
2im 2 -^-- 
E L 



2j3D d (2B 2 



+2m,n 



- d f QB 2 

uj— In 

ds \ E 



nw d / W£\ 

2 — dS{ ln ^l) 



(13) 



(14) 



We multiply both sides of this equation by $* M (where 
the star indicates the complex conjugate) and integrate 
over s between boundaries s m i n and s max . After an 
integration by parts we get a quadratic form: 



ds r 2 ( S )$; / ( S )i3f( S ) 



2Bl 
L E 



d * 

OS 



3. Corotation and Amplification 

3.1. A variational principle 

In this section we derive a variational principle which will 
allow us to analyze the effect of the corotation resonance, 
and also to discuss the boundary conditions to be used. As 
usual with variational formulations, this permits a more 
rigorous and general result than in previous discussions of 
both the resonance and the boundary conditions. In order 
to minimize the complexity, we present here only a simpli- 
fied form retaining only magnetic stresses; the full compu- 
tation, including pressure forces, is given in Appendix A. 
We start from Eqs. (0-^), which give: 



U = 



-2lUJ^ — <f> M 

E os 



d 



If Bn _ 

v= L\- 2W lta-s*» 



UmQ—<f> M 



(11) 
(12) 



2B 2 



EL 2 



i 2 {L + Ann') 



-2mn 



flW d 

uj ds 



hi 



WE 

"r 2 



2Bj 
L E 



M 



ds 



4> 



M 



(15) 



Let us consider the first term in the integral, in the 
left-hand side of this equation: using (||) we can write it 
as: 



dsr 2 B? (s) / ds'r'K m - Bf*(s') 



Using the property that: 

r 



i 

= -K„ 

x 



(x) 
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we see that this expression is hermitian (symmetric to the 
exchange of B° and B^* ) . 

Let us now consider the other terms in the left-hand 
side of Eq. (|l5|) : they have resonant denominators at L = 
(the Lindblad resonances, u> — ±/t) and at Co — (the coro- 
tation radius). We will admit without demonstration the 
common result (obtained by a Frobenius expansion) that 
the solutions are regular at Lindblad resonances (Narayan 
et al. , 1987), but not at corotation where the waves can 
exchange energy with the flow. This resonance has been 
discussed by Pannatoni (1983) in the context of galactic 
spiral waves, and by Papaloizou and Pringle (1985) and 
Narayan et al. (1987) in the context of the Papaloizou- 
Pringle instability. Its effect on the magnetically driven 
spiral density waves is the main object of this paper. 

The corotation resonance appears only in the last term 
in the integral, on the left-hand side of Eq. (15). As will 



In section 3.2 we discuss the main properties of the unper- 



be seen later, it contributes an imaginary term to the in- 
tegral and thus makes it non-hermitian. We will consider, 
as in the above-mentioned treatments, that this term is 
small (this can be obtained by reducing the derivative of 
Bq/WE), so that we can treat it as a perturbation. Then 
to lowest order, all the terms in the left-hand side are Her- 
mitian, and the right-hand side is a boundary term (con- 
taining the physics associated with boundary conditions, 
as will be discussed below). This means that Eq. flU) is 
self-adjoint, so that it forms a variational principle (Morse 
and Feschbach, 1953). Appendix A shows that, as could 
be expected, this property is still valid when we retain 
pressure forces. 

For readers not familiar with variational principles, we 
will only mention here their main advantage, namely that 
they allow a rigorous treatment of perturbation problems: 
let us note the variational form as: 

<r|£HIO = o 

where C = Cq + eC\ , e is a small parameter and £ is the 
vector formed from the components (velocities, potential, 
density) of the perturbation. We assume that we know 
(u>o, £o); an eigenvalue and eigenvector (the frequency and 
spatial dependence of the solution) of the unperturbed 
problem Cq: 

£o(^o)£o = 

and look for solutions (uiq + Slu, £o + &;) of the perturbed 
problem, to first-order in e. The hermiticity of Co gives 

(&\£oM\5S) = (6e\C (tj )\Z ) = 

so that 5lo can be obtained to first order without need to 
compute 5£, making the task relatively simple and more 
straightforward than the "direct" approach used e.g. by 
NGG. One gets simply: 



dui 



&)*w = -e<f5IA(wo)|fo> 



(16) 



turbed solutions. Since this neglects the corotation reso- 
nance, they are essentially the properties of magnetically- 
driven spiral waves, which had been obtained in the 



shearing-sheet approximation in Paper I. In section 3.3 we 



will perturb them by the corotation resonance term. The 
boundary conditions, appearing in the right-hand side of 
Eq. ( |l5| ) , are discussed separately in Appendix B and im- 
plemented in the numerical solution presented in section 

3.2. Waves and modes 

In this sub-section we discuss waves and modes, as they 
can be obtained from a WKB approximation in the un- 
perturbed state (without the corotation resonance). The 
WKB approximation applies for spiral waves when the ra- 
dial wavenumber is much larger than the azimuthal one 
(giving tightly wound spirals). Then the first two terms, 
on the left-hand side of Eq. ( |l5| ) are larger than the other 
ones, which are proportional to m. We look for a local 
radial wavenumber k: 




d. 



M 



zfc(s)4> 



A I 



and use the relation between $ m and B 1 ? 



B 



D 



M 



z=0+ 



together with Eq. (0) which gives: 



d_ 

dz 



\k\ 



M 



A I 



where the minus sign comes from the constraint that <E>m 
decreases at z — > oo. Then from the first two terms in 
Eq. (n5h we get the dispersion relation: 



~2 = K 2 + 2£f N 



(17) 



where for completeness we have added the pressure term, 
derived from the full variational principle in Appendix A. 
This is the dispersion relation derived in paper I in the 
shearing sheet approximation. It is identical to the dis- 
persion relation of self-gravity driven spiral density waves, 
with the magnetic term acting as a negative self-gravity. 

From this dispersion relation one gets the basic prop- 
erties of these waves: one can have both leading (k < 0) 
and trailing (k > 0) waves, resulting in spiral structures 
wound in opposite directions. Their group velocities, on 
both sides of corotation, are respectively toward and away 
from the corotation radius (where u> vanishes). 

The waves propagate only where Eq. (17) has a posi- 
tive root for |fc|, i.e. for Co 2 — k 2 > (beyond the Lindblad 
resonances). In the corotation region they are evanescent. 
Thus a leading wave traveling from the center toward coro- 
tation will be reflected at the Inner Lindblad Resonance 
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(ILR) as a trailing wave traveling back toward the cen- 
ter. If the boundary condition at or near the center is, as 
usually considered, such that this trailing wave is reflected 
again as a leading wave, one can establish a cavity between 
the center and the ILR; for a discrete set of frequencies, 
given by an integral phase condition, the resulting lead- 
ing wave has the same phase as the initial one, so that a 
standing pattern can be established. These patterns are 
known as normal modes of the system. 

Another property of the waves is that they have a neg- 
ative energy inside corotation, because they rotate slower 
than the gas, and positive energy outside. If the waves can 
tunnel through the forbidden band at corotation, they are 
amplified: let us explain this by considering a leading wave 
of energy -1 traveling from the center to the ILR, resulting 
by tunnel effect in the emission beyond the Outer Lind- 
blad Resonance (OLR) of a trailing wave of energy +e, 
unspecified but positive. Then the trailing wave reflected 
to the center must have energy — 1 — e, i.e. a higher am- 
plitude than the incoming wave: as it is reflected at the 
ILR, the wave is amplified. 

In the case of the Papaloizou-Pringle instability, where 
the only force comes from the pressure (which acts only 
locally), e is exponentially small: this is natural since the 
waves are evanescent in the forbidden band, as usual with 
tunnel effects, so that the energy of the transmitted wave 
is proportional to the exponential of a large negative quan- 
tity, of the order of —rfl/cs- However we have shown in 
PTS and paper I that the long-range nature of the forces 
change this result in the case of self-gravity driven or mag- 
netically driven spirals. This appears mathematically, in 
the shearing sheet model, as the contribution of a branch 
cut in the complex-fc space, related to the spatial behav- 
ior of the potentials. The reason is in fact quite simple: 
consider a density perturbation in the cavity between the 
center and the ILR. The tunnel effect allows it to act be- 
yond the forbidden band, in the same manner as in the 
Papaloizou-Pringle case. However, this density perturba- 
tion creates a potential which is not evanescent in the for- 
bidden band! As usual for an m-polar potential, it varies 
as r~ m , so that it decreases only algebraically across the 
corotation region. This radial dependence corresponds to 
the fact that the branch cut found in fc-space starts from 
singular points at k = im, or fc 2 
geometry of the shearing sheet. 



iky in the cartesian 



The result of this anomalous tunneling is that, for self- 
gravity driven spirals, e can be very large, rather than 
exponentially small as in the Papaloizou-Pringle case. 
In the magnetically driven case, we found in Paper I 
that maximal amplification is shifted to smaller azimuthal 
wavenumber, but remains small. We show in the next sub- 
section that the corotation resonance can make it much 
larger. 



Im(s) 



Integration contour 



— x — 



Re(s) 



-Pole at Q. (s) = co / m 



Fig. 1. The integration path for the integral in Eq. (IX Sf) . The contour 
lies on the real axis (giving a real Principal Part), except tor a clockwise 
semi-circle around the corotation pole. 



3.3. Resonant effect at corotation 

We can now turn to the perturbative treatment of the 
corotation resonance. For simplicity we will neglect bound- 
ary terms (the right-hand side of Eq. (|l5|)), treated sepa- 
rately in Appendix B. We thus rewrite Eq. (£L5[) as: 



ds r 2 ( S )^ f ( S )i?f( S ) 



2Bl 
L S 



os 



2B, 



o 



m r 8mil 2 W B„ 
ds 



EL 2 

,2 d 



$m| \ 0(m 
WE 



EL 2 



(18) 



where 0{m) contains terms we neglected in the WKB 
approximation; we will also neglect them here as in the 
"cavity" where the waves propagate (so that <E> m is large) 
these terms are small. To lowest order {i.e. neglecting the 
right-hand side ) this equation gives a discrete spectrum 
of normal modes, discussed in the previous sub-section, 
which are neither amplified nor damped (since the opera- 
tor is hermitian, oj is real). 

To the next order we perturb these solutions by the 
resonant term in the right-hand side , using Eq. ([l6|). We 
are only interested in the imaginary part of the integral in 
the right-hand side , since it will give us the growth rate 
resulting from the resonance. This is done in the classi- 
cal way used to study Landau damping in plasmas, or 
similar fluid problems (Lin, 1955, Narayan et al. , 1987), 
and for completeness we will sketch here its justification. 
One must remember that, since the modes we are look- 
ing for are exponentially growing or decreasing, the use of 
a frequency uj cannot be justified by a Fourier but by a 
Laplace transform; its inversion is obtained by an integral 
on a path in the complex-w plane, which must lie above 
{i.e. at larger imaginary part than) any singularity. In 
Eq. ( p"8| ) the singularity is at s such that oj = mfl{s), and 
we must compute the integral in the right-hand side as if 
u had a positive real part. One finds easily that this means 
that the pole s cor lies below the real-s axis, so that the 
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integration can be done on the contour shown in Fig. |l|: 
away from the pole the path stays on the real-s axis, so 
that its contribution is the principal part of the integral: 
this is essentially real, so that it only contributes to a small 
change of the real part of u. 

The rest of the integration contour is a half-circle 
around the pole, described clockwise because the pole 
(when the inverse Laplace transform is properly defined, 
i.e. when u> has a positive imaginary part) is below the 
real axis: its contribution to the integral is thus, from 
Cauchy's theorem, — in times the residue of the integrand 
at the pole. This gives: 



5u> 



ds 2 



0_ 

ds 



du/ds 



8m- 



n 2 w Bl 



d_ 

1$ 



or: 



du> 



ds 2^° 
S 



as 



-2uj 



d , (WY, 

7T ln 
as 



Uo 2 



(19) 



where the square bracket in the right-hand side is evalu- 
ated at the pole. The main contribution to the integral on 
the left-hand side comes form the inner cavity, at Q < 0. 
This integral is thus real and positive. On the right-hand 
side , fi' is usually negative, so that we get the sign of the 
growth rate 7 = Im(u;) resulting from the resonance: 



sign(7) = sign 



9-s ln {-Bj) 



(20) 



Thus outward gradients of W and £ are stabilizing, while 
an outward gradient of Bo is destabilizing. 

The full variational principle given in Appendix A con- 
tains a second resonant term, associated with pressure 
stresses. This term is proportional to the derivative of the 
specific vorticity (sometimes called vortensity in this con- 
text): 



sign( 7 ) 



sign 



J rcssurc 



d_ 

ds 



In 



(21) 



There is also a cross term, proportional to hB. Thus when 
magnetic stresses are absent we recover the usual result 
(Papaloizou and Pringle, 1985; Narayan et al. , 1987) on 
the effect of the corotation resonance. Note that in this 
case an outward gradient of £ is destabilizing. The result 
in the general case is a combination of these two contri- 
butions. 

However another result of the general variational form 
is that the pressure contribution is proportional to |<r| 2 ; 
on the oth er ha nd, and in spite of the parallel mentioned 
in section (2.3) between a and B®, the magnetic contri- 



bution is rather proportional to |$m| j where $m is the 



integral of B® ■ As mentioned in section (3.2), a is evanes- 
cent in the forbidden band at corotation, making the effect 
of the resonance exponentially small in the pressure-driven 
case. But although B® is also evanescent, $m is not! 

In order to understand this, let us imagine that the 
perturbed field B^ is strictly zero in the forbidden band, 
so that $m has its source only in the inner cavity where 
waves can propagate. Then in the corotation region <J>m 
will have the classical behavior of an m-polar potential in 
vacuum, varying as r~ m . This allows it to retain a siz- 
able effect at the corotation radius, and thus to efficiently 
exchange energy with the Rossby vortex. 

It would be of course difficult to find a simple explana- 
tion to the stabilizing or destabilizing effect of the corota- 
tion resonance, in this case as in the Papaloizou-Pringle in- 
stability. The starting point would certainly be that, when 
the equilibrium field varies across the corotation radius, it 
results in a variation of the tangential force exerted on the 
disk by the magnetic potential $m, produced by fluctu- 
ations in the inner disk region. This variation (or, in the 
opposite direction, that due to a gradient in the disk in- 
ertia) in turn results in differential azimuthal motions, at 
the source of the Rossby vortex. 

We will not try to get, from the variational principle, 
an estimate of the resulting growth rate. Rough estimates 
show that it should be of the order of va/t ~ /3 _1//2 fi/i/r. 
We will rather turn now to numerical solution of the equa- 
tions to confirm our analysis and exhibit explicit solutions. 

4. Numerical study 

4-.1. The eigenvalue problem 

We solve numerically the set of Eqs. ([?] - |o|), for the four 
unknowns U, V, h and B® , by turning it into an eigen- 
value problem. To do this we write the equations on a 
staggered grid, with velocities measured halfway between 
density grid points. 

The magnetic potential appearing in the right-hand 
side of the equations is expressed through Eq. (^) . Special 
care is needed for the Poisson Kernel, which is divergent 
at x = 1. We use a smoothing, replacing Eq. (||) by: 



K m (x) 



d$ cos md 



I" Jo (x 2 - 2ircosi9 + 1 + e 2 ) 1 /2 



(22) 



where e is of the order of cs/rtt, representing the effect of 
finite disk thickness; we have checked that varying e does 
not affect the results presented here. Then, for a grid of 
tig points separated by ds we define a vector £ of length 
4nc such that 



U ( s = (i + -)ds 



£(4») 



£(4i + l) = V[s=(i+-)ds 



Tagger and Pellat: Accretion-Ejection Instability 







£(4* + 2) 
£(4* + 3) 



h(s — i ds) 
B? {s = i ds) 



with i = 1. . . uq. In this manner our equations can be set 
as an eigenvalue problem: 



(23) 



where A4 is a (Ancx-Anc) matrix. It would be band 
diagonal if we were dealing with a differential problem, i. e. 
local forces. But the magnetic potential (i.e. long range, 
non-local forces), which turns it into an integro-differential 
problem, makes M. a full complex matrix, for which we 
seek eigenvalues and eigenvectors. M. would be hermitian 
in the absence of the radiating boundary condition at large 
s, and of the corotation resonance. 

One usually considers that 1000x1000 is the maximal 
size of a matrix that can be inverted numerically without 
excessive accumulation of roundoff errors. This leaves us 
with a maximum uq — 250: in practice, using the complex 
integration path described below, we have found this limit 
largely sufficient to get very good convergence of eigenval- 
ues and eigenvectors, and we have not tried to challenge 
it. The determination of eigenvalues and eigenvectors is 
thus done with standard library packages and takes only 
a few minutes on a modern workstation. 

4-2. The complex contour in s 

The main limit to the precision comes from the vicinity of 
the corotation radius, where the solution varies rapidly- 
the more as we consider very weakly amplified solutions, 
so that the pole lies very close to the real-s axis. We solve 
this difficulty, and at the same time the outer boundary 
condition, by solving on an oblique axis in the complex-s 
plane. 

This means that, taking s m i n = without loss of gen- 
erality (since s is logarithmic), we choose s max complex 
with a positive imaginary part. This serves at the same 
time three purposes: 

1. It allows us to move the integration axis far from the 
pole, in order to get a smoother behavior. Given the 
limited number of grid points we can use, this strongly 
improves the accuracy of the eigenvalues. On the other 
hand it allows us to exhibit the eigenvector only on this 
complex axis, so that its components (the velocities 
etc.) are not the physical solutions one would obtain 
by solving on the real axis. However we maintain the 
obliquity of the integration axis weak enough to limit 
this problem. In selected cases, where the growth rate 
is strong enough, we solve on the real-s axis and find, 
as could be expected, that the solution is significantly 
modified only at large s (note that our solution on the 
complex axis is correct, but is not the one physically 
observable in real space). 



Fig. 2. A typical profile of the magnetic field (solid) and its logarith- 
mic derivative (dots). We have a strong but localized gradient, resulting 
in a limited total change of Bq. A similar profile is used for the surface 
density of the disk 



Given the constraints on the inverse Laplace trans- 
form, 



3.1, solution along this 



discussed in section 
oblique axis allows us to correctly describe damped so- 
lutions by analytic continuation, as long as the imagi- 
nary part of Smax is sufficient that the integration con- 
tour lies on the correct side. We will see below that this 
appears in a very obvious manner in the eigenvalue 
spectra. 

The boundary condition we apply at s max is the ra- 
diation condition, also called the outgoing wave con- 
dition. This condition states simply that we do not 
want the system to receive information from outside, 
so that the only perturbation allowed at s ma x is the 
wave traveling outward. In our case, far enough be- 
yond corotation that the WKB condition is valid, this 
translates into accepting only a trailing wave, with a 
loca l radial wave number k > 0. We show in Appendix 
|B.l| that this same obliquity of the integration axis 
(with Im(s > 0)) provides a simple and efficient way 
of implementing this boundary condition. 



4-3. The profiles 

We wish to describe an equilibrium with gradients of den- 
sity and vertical magnetic fields. On the other hand we 
do not wish too strong changes of these quantities, which 
would cause large differences in the radial wavelengths and 
thus create numerical problems. Thus we use the analyt- 
ical result, showing that amplification depends only on 
the local gradients at corotation. This allows us to take 
profiles with strong but localized gradients, giving limited 
total variations across the integration domain: 



E(s) = S [1. - d s taxih (p E (s- s ))] 
B (s) = B 0int [1. - d B tanh(p B (s - s ))] 



(24) 
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Using or 4 ~ .2 gives reasonable profiles for our pur- 
pose. This does not claim to represent realistic profiles, but 
allows us an optimal study of the physics involved. Fig. (||) 
shows a typical profile of the magnetic field and its loga- 
rithmic derivative. When needed we use trial and errors, 
on low-precision runs with a limited number of points, to 
change sq so that the corotation radius lies close to the 
maximal gradient. 



4- 4- Numerical results 

We measure the strength of the magnetic field by the usual 
parameter giving the ratio of thermal to magnetic pres- 
sure: 



(3 



8np 



In B 



1 



taken at the inner disk radius. At low field > 1) we 
expect the magnetic component of the perturbation to 
play only a weak role. Furthermore, this is the regime 
where the magneto-rotational instability is unstable and 
should dominate anomalous viscosity. Nevertheless we will 
present numerical results at (3 — 50, in order to check the 
general physics presented here. In contrast we will show 
a second set of results for = 1, where our instability is 
obtained while the magneto-rotational instability (which 
does not appear in our analysis since we consider infinitely 
thin disks) should be stable. 

We use the result of Eq. ([l9]) to limit the parameters 
we vary: we keep £ constant and introduce only a radial 
variation of Bq . We use trial and error to have a maximum 
gradient very near the corotation radius, with 

d_ 

ds 

so that, from Eq. (|l9|), the modes should be weakly unsta- 
ble if the magnetic contribution dominates, and (from the 
additional terms included in Eq. ( |A.3| ) ) damped if pres- 
sure dominates. We choose a sound speed cs — .lr rt , 
where the subscript notes values at the inner radius 
(.s = 0). Figs. (§-^) show the spectrum of eigenvalues 
we obtain for /3 = 50 and m = 4. Frequencies are nor- 
malized with f^o = 1., so that eigenvalues whose real part 
lie between and m have a corotation at s > 0. We 
distinguish three components in the spectrum: 

1. a set of values aligned diagonally, as seen in Fig. (|J). 
From their behavior as the number n s of grid points is 
varied, we find these values to form a continuous spec- 
trum as n s — y oo, associated with the discretization of 
the problem, and thus not physical. 

2. A second component, labeled (a) in Fig. (^). These 
also form a continuous spectrum, and correspond to 
unphysical solutions with corotation very near grid 
points: thus when we solve along a real-s axis, this 
component lies on the real-w axis. This permits us to 
directly check on the plots that all roots above this 




Fig. 3. The full spectrum of eigenvalues in a typical low-field case, 
with j3 — 50. and m — 4. Integration is performed on an oblique complex 
axis s — (0., 0.) to (1.5, .1). Most of the eigenvalues are associated to 
the discretization of the system. 




Fig. 4. A zoom on the results of Fig. Most noticeable are (a) a dense 
component corresponding to a continuous spectrum with corotation on 
the integration axis (b) another component corresponding to the physical 
discrete spectrum, close to Im(oj)— 0. 



set, even when (as here) they are weakly damped, lie 
on the proper side of the integration axis, as discussed 
in section |f.2| . 

A third component, labeled (b) in Fig. (^), corresponds 
to a discrete spectrum. In order of decreasing lu, these 
eigenvalues correspond to the classical set of modes 
with n = 0, I, 2, . . . nodes in their radial structure, be- 
tween the inner radius and corotation. Due to the lim- 
ited numerical precision, the frequency of only the fun- 
damental solution, n — 0, is reliable, although higher- 
order ones do show the increasing number of nodes, 
verifying our identification of this discrete spectrum. 
The n = eigenvalue converges to a very good accu- 
racy when n s is varied above typically 200. 
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Fig. 5. Real (full line) and imaginary (dashed) parts of the n — tti- q 

eigenvalue for jS = 50. All the modes are damped by the corotation * 1 S" 8 - The frequency and 



resonance. 



m = 2 mode. 



growth rate as a function of 0, for the 




Fig. 6. The spectrum for ft — 1, m — 3. The fundamental mode 
has a significant positive real part, i.e. is amplified by the corotation 



resonance. 




Fig. 7. Real (full line) and imaginary (dashed) parts of the n — 
eigenvalue for j3 — 1. The modes have become unstable by the corotation 



resonance. 
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Fig. 9. The eigenvector: Real (solid) and imaginary (dots) parts of U, 
Real (dashes) and imaginary (dash-dots) parts of V for /3 — 1, m — 3. 
Corotation is at s — .28. 



Fig. ([s]) shows the variations of the pattern frequency, 
uj/m, and growth rate, 7 — Im(w), as a function of m. 
At this high value of (3, the modes are all weakly damped 
by the corotation resonance, according to Eq. (fH|). 

On the other hand, Fig. (^|) shows the spectrum for 
j3 = 1., m = 3. The fundamental mode has acquired a 
significant positive imaginary part, i.e. becomes unstable 
by the corotation resonance. Fig. (^) shows that low-m 
modes become unstable. The growth rate is not strong 
{e.g. for m — 2, 7 is .06 times the pattern frequency), 
but one has to remember that we are dealing with normal 
modes, i.e. long-lived, exponentially amplified structures, 
so that this may result in strong amplitudes after a finite 
time. 

We show on Fig. (|^) the mode frequency and growth 
rate for m = 2, as a function of (3. The frequency decreases 
as (3 decreases, while the growth rate has a maximum for 
P close to 1. This can be understood from the dispersion 
relation, Eq. ([l7]): at low f3 the wavelength becomes larger, 
widening the "cavity" and thus pushing to larger r the 
corotation radius ; on the other hand the turning point 
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Fig. 10. The eigenvector: Real (solid) and imaginary (dots) parts 
of a. Real (dashes) and imaginary (dash-dots) parts of <&m for — 1, 
m = 3. Corotation is at s — .28. 



(the radius where the radial wavenumber k s = 0) is in 
fact not exactly at the Lindblad resonance, but where 



oj 2 



k = h 



(25) 



so that the forbidden band around corotation widens as 
P decreases. Thus at low (3 the amplification mechanism 
becomes less efficient. 

Finally we show on Figs. (| - 0) the eigenvector for a 
typical mode, showing in particular the cavity inside coro- 
tation and outgoing wave beyond it, decreasing at large s 
because of the obliquity of the integration axis. 

5. Discussion 

We have presented an instability of magnetized accretion 
disks, which has the important potential of relating accre- 
tion (by a spiral density wave in the disk) and ejection, 
through Alfven waves propagating to the corona of the 
disk the energy and angular momentum extracted by the 
spiral. The instability has the following characteristics: 

— It forms normal modes, i.e. standing, exponentially 
growing, wave patterns. 

— Their growth rate is rather small, typically of the order 
of hfl/r where h is the disk thickness (we could have 
shown stronger growth rates by increasing the gradi- 
ents in our numerical examples). But their exponential 
behavior still allows them to reach high amplitudes in 
a finite time. 

— The modes are localized in the innermost part of the 
disk: since the sound and Alfven velocities are small 
compared to the rotation velocity, the wavelength in 
the "cavity", between the inner radius and the ILR, 
is of the order of a few times the disk thickness (we 
will not discuss here the exact relationship, which is 
not important). Thus the frequency of the fundamental 
mode (with no node in the cavity) is essentially fixed, 



at a given m, by the existence of an ILR close to the 
inner disk radius: lu — mfl — —k at r — Ti close to the 
inner radius tq, with il = k in a keplerian disk, gives 
uj = mQ(r c ) = (m— l)O(ri), where r c is the corotation 
radius, and thus 

2/3 



i - 1 

Thus these modes cannot be considered as candidates 
to explain accretion in the bulk of the disk, but only to 
energize jets close to the inner radius, where, in usual 
viscous models, half of the energy and momentum of 
accretion are deposited in a boundary layer. Actually 
the size of the "cavity" is such that it might be con- 
sidered as a form of boundary layer. 

— The instability occurs in strongly magnetized disks, 
with (3 ^ 1 (where the magneto-rotational instability 
is stable) 

— It appears when the quantity £IY,/Bq increases radi- 
ally. This requires either an inverted density gradient, 
or a strong magnetic field gradient. This means that 
disk models must consider the radial transport of mag- 
netic flux, as well as that of the accreted matter: the 
horizontal part of the magnetic flux can be expelled 
from the disk vertically by buoyancy (e.g. the Parker 
instability, Foglizzo and Tagger, 1994, 1995), but the 
vertical part can only accumulate in the central region 
of the disk, if the accretion turbulence transports it to- 
gether with the matter. One may note as an example 
that at the center of our galaxy, on a scale of ~ 100 pc, 
the magnetic field seems to be vertical and of the order 
of a few milligauss, as compared with a horizontal field 
of a few /j,gauss elsewhere in the galactic disk. 

This threshold in radial gradients also implies that 
this instability is very likely to result in a bursty behavior 
of the accretion/ejection process, associated with the 
processing of both gas and magnetic flux in the disk. 

We also wish to mention the possible connection be- 
tween this instability and the QPOs observed in X-ray 
binaries: since the mechanism favors low m values, i.e. 
spirals with a small number of arms, it is very likely to 
appear as very coherent, monochromatic structures, in the 
manner of galactic spirals. Their corotation radius would 
be of the order of a few times the inner disk radius. These 
are precisely characteristics of the low frequency QPOs, 
with a frequency in the range of .1 (in black hole bina- 
ries) to a few tens (in neutron star binaries) of Hertz (see 
e.g. Wijnands and Van der Klis, 1999), i.e. a frequency 
typical of keplerian rotation not far from the inner disk 
boundary (Markwardt et al. , 1999 ; Psaltis et al. , 1999). 
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Appendix A: The full variational principle 

In this append ix w e perform essentially the same computations 
as in section ( ]3.l[ ), but retaining pressure stresses to give a 
complete variational principle. 
Thus Eq. @ becomes: 
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In the same manner, starting from Eq. 
similar equation for a: 



(A.l) 
we obtain a 



2 9 (2B 9 , \ 9 (Ed 

ra = ~dS ^ M J _ 9s" lL9 S - 

+^ < ^>A/<{™ 2 (L + 4QQ') 



c s h 



+2mfl 



9s 



cD 9s 



+ -^c|/i<j m 2 (L + 4fifi') 



L 2 



- 9 , „W 9 , (W 

lo — ln(nE) +2 —In 

9s cj 9s 



(I) 



(A.2) 



We multiply Eq. (|aJ) by c§ft*, and Eq. by 2$^, 

sum the results and integrate over s: this gives, after integra- 
tions by parts, the full variational principle similar to Eq. (Ila): 
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As in Eq. ([lE]) the left-hand side is hermitian, so 
is again a variational principle. 

Appendix B: Boundary conditions 

B.l. Outer boundary: outgoing wave 



(A.3) 
that this 



As explained in section 4.2, we require that at the outer bound- 



ary the solution reduces to an outgoing wave. The complex 
integration axis allows us an easy implementation of this con- 
dition. We start by assuming that s max is large enough that 
the WKB approximation can be used. 

Then let us consider the two WKB solutions: 



ds'k(s' 



f± - exp ±i 



where the plus and minus sign correspond respectively to the 
trailing (outgoing) and leading (incoming) waves. If we inte- 
grate on an oblique axis in the complex-s plane, with Im(s) 
positive, one sees that the trailing solution decreases exponen- 
tially, while the leading one grows exponentially, as s increases. 

Let us now assume that we have at s max an incorrect 
boundary condition. In general, this means that we have lead- 
ing and trailing waves of comparable amplitude. But as one 
goes inward from s max , the trailing solution grows exponen- 
tially as s decreases while the leading one decreases. Thus at 
some distance inward from s ma x the correct (trailing) solution 
will dominate. This method of dealing with the boundary con- 
dition has been introduced in disk physics by Drury (1980). It 
allows us here to take any simple (but a priori false) boundary 
condition at s ma x , and only check a posteriori that the correct 
behavior results; if it is not the case we just need to increase 
lm(s max ). In practice we use a reflecting condition at s max and 
have little difficulty with the outer boundary condition . 

In fact we have another effect helping us with the outer 
boundary condition: as seen from the WKB dispersion rela- 
tion, Eq. (^), the wavelength becomes small at large s. On 
the other hand our numerical scheme for the discretization of 
the equations is such that, when the wavelength becomes com- 
parable with the grid spacing, waves suffer a heavy numerical 
damping. Then if we have a solution resulting in an outgoing 
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wave, this wave is damped at large s and very little of its energy 
can be reflected. Thus taking Re(s max ) large enough is also a 
good way of implementing the outer boundary condition. 

Let us now turn to the contribution of the outer boundary 
to the variational principle, Eq. (|l5|). It is easily obtained in 
the WKB approximation, 



d_ 

ds 



with k given by the dispersion relation, Eq. (J17|) . Let us treat 
this boundary term by perturbations, as we do for the reso- 
nance term: we see that, since it is imaginary, its contribution 
will be to the growth rate of the mode. In fact one easily checks 
that with k > (for an outgoi ng w ave) it is always destabi- 
lizing, as mentioned in section ( |3.2|) . This is the amplification 
mechanism discussed, in the shearing sheet approximation, in 
paper I. The resulting growth rate is always quite small. 



B.2. Inner boundary: reflection 

Any condition at the inner edge of the disk, which does not 
change the energy of the wave, will reflect it with an efficiency 
of 1. This is the classical picture leading to normal modes. 
Various conditions satisfy this criterion, and we discuss them 
here with the help of the variational principle. For simplicity 
we will discuss here only the pressureless variational principle 
of Eq. ( |l5| ) , but the generalization to the complete problem is 
straightforward. 

From Eq. we see that by making either 



or 

d_ 



o 



§M = o 



(B.l) 



(B.2) 



the boundary term vanishes, so that the mode energy is not 
changed. Either of these conditions is thus an acceptable re- 
flecting boundary condition. 

On the other hand it is more usual to assume a rigid bound- 
ary and require that 



U = 

at the inner 
terms, give: 



(B.3) 



edge. In this case Eqs. (fjjj^) , forgetting the pressure 



d 



R 2 

4im-|$^^$ M = 2iu*Q,\V\* 

Zj^ OS 



(B.4) 



making the boundary term real. Thus if we start from a 
marg inal mode (u> real) with boundary conditions (B.l) or 
( |B.2| ), and change it to U = 0, the resulting change to lo is 
real: this means that condition (B.3) also gives a perfect re- 
flection of the wave, but with a different phase; the resulting 
change in the integral phase condition modifies the discrete set 
of mode frequencies, but does not introduce growth or damp- 
ing. We note that the condition V — use d by NGG would 
give a similar result. We use condition (B.3) in the numerical 
solutions. 



Appendix C: Emission of Alfven waves in the 
WKB limit 

In this appendix we describe in the simplest manner the emis- 
sion of Alfven waves, carrying some of the energy and angular 
momentum of the density waves in the disk to its corona if the 
latter has a small but non-vanishing density (thus relaxing the 
hypothesis, made throughout the paper, that the disk lies in 
vacuum) . We do this in a simple WKB limit in the region where 
density waves propagate radially, ie far from the corotation re- 
gion. A full computation, without the WKB approximation, 
would be necessary to describe the emission of energy and mo- 
mentum from the vortex at corotation. We defer it to a future 
publication. 

In order to describe the vertical propagation of the Alfven 
waves, we now consider the z dendence of the variables. We 
start from the Euler equations and derive equations for two 
new quantities: 

D(s,z) = -^-U + imV = r 2 Vi • V± 

OS 

R(s,z) = ^-V - imU = r 2 V ± x V a_ 

OS 



where the subscript _L refers to the components of a vector in 
the plane of the disk. We get: 



iu>D - 2QR - 20,' (v - imu) 
-iuR + WD + W'u 



= r 2 V ■ (F/po) 
r 2 V x {F/p ) 



(C.l) 
(C.2) 



where po(s,z) is the equilibrium density, and F is the per- 
turbed force acting on the fluid. Here for simplicity we will 
again neglect pressure forces, unimportant for Alfven waves. 

We now use the WKB approximation, by defining a local 
radial wavenumber 

~ ik 

OS 

with k 2> m. The dispersion relation, eq. (p^), tells us that this 
is verified far from corotation. Consistently neglecting all radial 
gradients of equilibrium quantities, and using the induction 
equation for the perturbed magnetic field, we get after some 
algebra: 



iuD - 20R 
-iujR + WD 



2 2 

no r £ 



dz- 



' dz 2 
■R 



D 



(C.3) 
(C.4) 



We now use the same procedure as in Tagger et al. (1992): 
let us assume that at large z (z S> h, where h is the disk scale 
height) the density goes to a small but non-vanishing density 
Poo mimicking a di sk corona. There va becomes very large, so 
that equation (C.3) gives to lowest order: 



D(z) 



-\kz\ 



where the minus sign comes from the condition that D de- 
creases at infinity. Thus we get the same behavior derived for 
$m in the cas e wh ere the disk lies in vacuum. On the other 



hand equation(C4) gives: 
R(z) ~ e lk * z 
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where 



~> Va = ~a 



and causality requires that k z has the same sign as Co, so that 
this describes an Alfven wave e mitted fr om the disk. 

We now multiply equations ( |C.3| - |C.^ ) by p(z) and integrate 
over z. We get: 



iuD - 2QR 
-iCjR + WD 



__A 

iui rh 



q D 



(C.5) 
(C.6) 



where the bar means an ave rage over the disk height, and the 
right-hand side of equation (C.E) is an estimate of the integral, 
which allows us to recover in the limit poo —* the disper- 
sion relation, equation jl7j). From these equations we get the 
dispersion relation when the coronal density is small: 



J*l. 
rh 



L0 2 h A 



(C.7) 



Using the value of k z , and solving by perturbation, we find 
that the emission of Alfven waves to the corona results in a 
damping of the density wave, i. e. uj gets a negative imaginary 
part: 



"7 'Alfven 



1 



2w 2 / 3 1 /2 



K 2 n f 



1/2 



(C.8) 



where po is the density at the midplane of the disk. From this 
result we can get the following conclusions: 

— The damping is proportional to the square root of the den- 
sity at infinity, a result typical of magnetic breaking. 

— It applies through R, the torsional component of the per- 
turbation: as expected, R involves a torsional (in contrast 
with D, giving a compressional) motion of the footpoints 
of the field lines in the disk. This motion propagates to the 
corona as an Alfven wave. It is noteworthy that R is sin- 
gular at corotation, because o f th e vorticity gradient con- 
tribution (W'u) in equation (C.2), and indeed represents 
the Rossby vortex there - although this is not described in 
the WKB approximation used in this appendix. Thus we 
can expect a strong Alfven wave emitted from thee Rossby 
vortex at corotation. 



This may be even stronger, since equation (C.J) shows a 



vanishing denominator at corotation. This must of course 
be taken only as an indication, since the WKB approxima- 
tion vanishes in this region. 
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The full computation of this effect would need to relax the 
WKB approximation. Its consequences on the formation of a 
jet requires non-linear study of the deposition of Alfven wave 
energy in the corona. These extensions of our result will be the 
object of future work. 



References 

Balbus, S.A., and Hawley, J.F., 1991, ApJ 376, 214 
Binney, J. and Tremaine, S., 1987, Galactic Dynamics, Prince- 
ton University Press 



